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Abstract 

A new way of diagonalizing the Jaynes-Cummings Hamiltonian is proposed, which allows 
the definition of annihilation operators and coherent states for this model. Mean values and 
dispersions over these states are computed and interpreted. 


1 Introduction 

The Jaynes-Cummings (J.C.) model |l|, which is extensively used in Quantum Optics, describes, in 
its simplest version, the interaction of a cavity mode with a two-level system. In the rotating-wave 
approximation, it may be described by the Hamiltonian [1,2] 

H iC = v (a j a + ^ a 0 + y<r 3 + /c(aV_ + aa + ) , (1) 


where {a,,a 2 ,a 3 } are Pauli matrices, cr 0 is the identity, er ± = o l ±ia 2 , and a ) and a are the photon 
creation and annihilation operators. Moreover, k is the coupling constant, to is the field mode 
frequency, and u 0 is the atomic- frequency. Let us also introduce the detuning A = u) — u> 0 . The 
exact solvability of this model is well-known. Working in the Fock space 
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the energy eigenstates take the form (for n — 0, 1,2, . . .) 
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In the expression of r(n) we have introduced the parameter 6 which will be important in the 
following. The corresponding energy eigenvalues are 

E~ = un + Kr(n), E+ = u(n + l) - Kr(n + l) . (7) 

The interest of this model, its solvability and its applications have long been discussed. More 
precisely, dynamical properties have been obtained through the use of states which are initially 
harmonic oscillator coherent states [3], but that evolve according to the J.C. Hamiltonian (2, 4]. 
Here, we construct new coherent states which correspond to eigenstates of an annihilation operator 
for H JC . To do that, we have to find first such an annihilation operator through the diagonalization 
of the Jaynes-Cummings Hamiltonian (1). Second, we use the theoretical approach, based on the 
direct product of the W6yl-Heisenberg group with SIJ(2), to evaluate those coherent states. Finally, 
we exhibit some of their properties. More details can be found in [5]. 


2 Annihilation operators and coherent states for H JC 

i 

The diagonalization of H iC is performed by the unitary operator O, so that 


H d = O'HkO = 


uj(N + 1) — kt(N + 1) 
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0 

uN + nr(N ) 


( 8 ) 


where N = (Aa, and the definition of r(N) is given in Eq. (6). This operator O has the form 


o = 


1 


R(N+ 1) L2 


+ Kr(N + 1) 


K 


— Cl' 


K 


1 


R(N + 1) 


a 


R{N + 1) 


R(N) [2 


A 

— + k r(N) 


( 9 ) 


Clearly, an annihilation operator for H d is given by A d = acr 0 . Since the states depend also 
on the spin index ±, we introduce the spinorial annihilation mid creation operators: 



We then obtain candidates to be annihilation and creation operators for H }C as: 


A = OA d O\ A j = OAlO f , E ± = OE ±d O f . (11) 

For the determination of coherent states, the situation is particularly simple when we work 
with H d in (8), since the energy eigenstates are the Fock-space basis vectors (2). The group 
theoretical approach to coherent states leads us to define 


\ z ,P) d =T(zJ]) d \ 0,-), z,pe C, (12) 

in terms of the unitary representation of the direct product of the Weyl-Heisenberg group with 

SU {2) 

T(z, 0) d = exp[^J - zA d + /iE +d - [m_ d ]. (13) 
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The coherent states for H ]C are then given by 

\z,ft) = O\z,0) d = T(z,(3)\EZ) = OT(z,ft) d O'\E»). 
From the harmonic oscillator coherent states we easily get 
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cos(0/2) | z—) d + e l<i> sin(0/2) |z+) d , (15) 


with ft = (0/ 2)e^, and | z) the normalized state \z) = e~ ^ The s<: ' ate \ z iP)d 

is a linear combination of the “fundamental coherent states” |z+) d and | z—) d , which are both 
eigenstates of A d (but npt of £_ d ). Similar fundamental coherent states are defined by using such 
a decomposition of \z,ft) in (14). 

To be complete, we give the time evolution of such states; the one of the general state (lo) will 
be obtained easily by linear combination. We start with the diagonal case, for which the evolution 
operator is 

/ gtt|w(JV+l)-isr(JV+l)l 0 \ 
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and we compute | z,±,t) = OU d (t)\z,±) d = U JC (t)0\z,±) d . More precisely, we get 
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These are similar to the states obtained by the evolution of the harmonic oscillator coherent states, 
except for the supplementary oscillation for each n in the sum. This implies that the coherent state 
does not evolve in time to another coherent state, unlike to the case of the harmonic oscillator. 
Moreover, our states (17) are different from those considered in other approaches [2, 4]. Indeed, 
all these authors deal with the states \z,±) d (or a mixture of them) and their evolution is 

U ]C (t)\z,±) d = e- itH “\z,±) d . (1«) 


3 Relevant physical quantities 


Let us recall the introduction of the parameter 6 in the expression of r(n) in (6). It will be used 
as a variable in the following. It contains both the detuning A and the coupling parameter k 
and leads to the exact resonance case when 8 = 0 or to the weak coupling limit when 8 — ► oo. 
The parameter x = \z\ 2 is also introduced, and will be proved to be a good approximation of the 
number of photons. We will deal with the function 


G(8,x) = e~ x 




(19) 


Its asymptotic behaviour is G(8,x) ~ y/x, which is independent of 8. The calculations will be 
done explicitly over the states (17), and for a general state | z,ft,t). We can use (15) to write 


383 



the mean value and dispersion of an operator X. Indeed, defining (X) ± = {z,t.,±\X\z,t,±) and 
(A r ) ±T = {z, t, ±\X\z, t, t), the mean value of X over a general coherent state is 


(X) = + i±|^ <*>- + ^(e* W+ _ + e-* W _ + ). 


( 20 ) 


In the case of having (X) + _ = {X)_ + =0, the square of the dispersion is simply 
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(AX) 2 = (X 2 ) - (X) 2 = ” — (AX) 2 + —-rA- ' (AX) 2 + ^XLX[(x) + - (X )_ 12 


( 21 ) 


Let us now compute the relevant physical magnitudes. The operator M — (a) a + 1/2)<t 0 + <x 3 /2 
corresponds to the total number of particles. It is a constant of motion and is invariant under the 
transformation by O. We then get 


(N) + =x+l, (N)_=x, (A M)l = (AH) 2 _ = a: . 


( 22 ) 


(Those are known results in connection with the susy harmonic oscillator.) 

The evaluation of the mean values of the number of photons N = a^a is less trivial, and gives 
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Comparing with the harmonic oscillator, where (N) = x, we have a correction due to the inter- 
action. Since k is usually small, a good approximation to the average number of photons is x. 
Indeed, the contribution of the terms containing the series is, in this case, approximately 1 /2. 

We can compute the mean values and dispersions of the energy in the fundamental states and 
study their behaviour with respect to both x and 8. Since (//jc) + _ = {//jc)_ + =0) the calculations 
over the general coherent states through Eqs. (20)-(21) do not give anything new with respect to 
the results for the fundamental states. Indeed, we can see from (21) that the dispersion attains 
his minimum over the pure states. The mean values are easily computed and take the simple form 


<//jc> + =w[(ar+l)-AG(«,x)], (H iC )_ = u[x + XG(8 - l,x) 
while the values of the dispersion are more complicated and present interesting features 


(24) 


(A H ic )\ = u 2 [a 2 (1 +8) + ( 1 + \ 2 )x + 2Ax(G(5,x) - G(8 +l,x))~ A 2 (G’(6,a;)) 2 ] , (25) 


(AHjc) 2 _ — w 2 A 2 S + (1 + A 2 )x — 2\x(G(8 — 1, or) — G(8,x)) - A 2 (G(8 — l,a:)) 2 . 


(26) 


We have introduced A = k/u. When a large number of photons is considered, we can use the 
asymptotic behaviour of G(8, x) to see that 


(A H JO ) ± 

(H , c ) ± 


(27) 


as in the harmonic oscillator case. 
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Fig. 1. Behaviour of (A//jc)+/w 2 as a function of 6 and x for A = 8. 

If we want to see how the dispersion evolves with respect to a variation of the characteristics 
of the system 8 and A, we analyse the form of (A H :c ) 2 + , since (A H JC ) 2 shows a similar qualitative 
behaviour. If we fix A, the typical behaviour of (A H ic ) 2 + is as in Fig. 1. It can be proved that for 
fixed values of 8 smaller than a certain 8 0 , (A H JC ) 2 + has a minimum for x f 0. 

The atomic inversion is the last quantity we will consider. Over a general coherent state, we 

have 
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where ip n (t) = <f> + 2 tnr{n + 1). If we take A = 0 and 0 = -(f) . 

behaviour which is similar to the one obtained by Narozhtiy ct dl. [2] (let us recall that their 
states are different from ours). It consists of Rabi oscillations, as shown in Fig. 2: 
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The derivative of this function with respect to t is essentially the value obtained in Ref. [2]. In 
Fig. 2, we show the graph of our (cr 3 ) for x = 20. It is similar to those obtained in many other 
papers, and that although the expression of (<r 3 ) is not exactly the same in all the cases. 

These results indicate that it is reasonable to analyze the coherent states associated to the 
.1.0. Hamiltonian in the way we are doing, but we do not know for the present if they could be 
interesting for the experiments. 
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x = 20 



nt 

Fig. 2. Collapses and revivals of the atomic inversion in a general coherent state for x = 20 and 6 = 0. 
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